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The density distributions of 10 Be and n Be nuclei obtained within the quantum Monte Carlo 
(QMC) model and the generator coordinate method (GCM) are used to calculate the microscopic 
optical potentials (OPs) and cross sections of elastic scattering of these nuclei on protons and 12 C 
at energies E < 100 MeV/nucleon. The real part of the OP is calculated using the folding model 
with the exchange terms included, while the imaginary part of the OP that reproduces the phase 
of scattering is obtained in the high-energy approximation (HEA). In this hybrid model of OP the 
free parameters are the depths of the real and imaginary parts obtained by fitting the experimental 
data. The well known energy dependence of the volume integrals is used as a physical constraint 
to resolve the ambiguities of the parameter values. The role of the spin-orbit potential and the 
surface contribution to the OP is studied for an adequate description of available experimental 
elastic scattering cross section data. Also, the cluster model, in which n Be consists of a n-halo and 
the 10 Be core, is adopted. Within the latter, the breakup cross sections of n Be nucleus on 9 Be, 

93 Nb, lsl Ta, and 238 U targets and momentum distributions of 10 Be fragments are calculated and 
compared with the existing experimental data. 
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I. INTRODUCTION 

The discovery of halo nuclei fjj has been related to the 
measured interaction cross sections of nuclei like 6,8 He, 
11 Li, Be isotopes with various target nuclei The 

evidences of the existence of an extended halo in neutron- 
rich nuclei are based on the observed unusually narrow 
momentum distribution of a core fragment and enhanced 
reaction cross section. Thu first example was the breakup 
of 11 Li at high energies [M by observing the large in¬ 
teraction reaction cross section .2} and the narrow mo¬ 
mentum distribution of 9 Li in the breakup of 11 Li, e.g., 
in the reaction n Li-|- 12 C at E = 800 MeV/nucleon Q. 
Here we should mention also the results of the experi¬ 
ments at lower energies (E = 60 MeV/nucleon) of scat¬ 
tering of n Li on 9 Be, 93 Nb and 181 Ta jll[ and of n Li on 
a wide range of nuclei from 9 Be to 238 U [l2|. As shown 
in Ref. (Ti|, not only scattering but also the breakup of 
11 Be in the collisions with the target nuclei 93 Nb, 181 Ta, 
and 238 U play a decisive role when studying the inter¬ 
nal cluster structure of 11 Be. Indeed, the narrow peak of 
the momentum distributions of the breakup fragments of 
such a neutron-rich nucleus reflects the very large exten¬ 
sion of its wave function, compared to that of the core 
nucleus 10 Be, and thus evidences the existence of the nu¬ 
clear halo [l4 - [20j . As was concluded in [jjj|, namely the 
longitudinal component of the momentum (taken along 
the beam or the ^-direction) provides the most accurate 
information on the intrinsic properties of the halo, being 
insensitive to details of the collision and the size of the 
target. In addition, recent measurements of the charge 
radii of 7 > 9 do,iig e pointed out that the average distance 


between the halo neutrons and the 10 Be dense core of the 
11 Be nucleus is around 7 fm (2l|. Thus, the halo neutron 
is about three times as far from the dense core as is the 
outermost proton because the core itself has a radius of 
only 2.5 fm. 

An important finding when investigating reactions 
with 10 Be and 11 Be nuclei, in particular the 10 Be+n 
breakup of n Be, is the effect of the deformed 10 Be core 
on the two-body cluster structure of n Be. In fact, in the 
n Be nucleus the inversion of the pi/ 2 and s i /2 orbitals 
predicted by Talmi and Unna [HJ and confirmed by Al- 
burger et al. [H] leads to a l/2 + ground state. Also, 
the probability of the El transition from this ground 
state to the 1/2“ first excited state of n Be located at 
320 keV excitation energy is the largest ever measured in 
light nuclei [U, [25| . The effects of the core deformation 
on the breakup of 11 Be on protons have been studied 
in several works. For example, in addition to a two- 
body cluster structure with an inert lo Be(0 + ) core and 
a valence neutron used in Ref. [26| in the continuum dis- 
cretised coupled-channels (CDCC) calculations of elas¬ 
tic and inelastic proton scattering on n Be, the authors 
have also discussed the necessity to account for contribu¬ 
tions from configurations involving excited states of the 
10 Be core to the 10 Be+n continuum of 11 Be. Crespo et 
al. [13] have found that the core excitation p+ lo Be(0/“) 
—> p+ 10 Be(2/") provides a significant contribution to the 
breakup cross section of n Be on the proton target at 63.7 
MeV/nucleon incident energy. 

In the earlier works (e.g., Ref. [28|) the elastic scat¬ 
tering cross sections of 10,11 Be on protons have been cal¬ 
culated using phenomenological OPs of given forms with 
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numerous fitting parameters of their real (ReOP) and 
imaginary (IrnOP) parts. However, in the further calcu¬ 
lations the more physically motivated microscopic folding 
models were applied (see, e.g., [HHIl]). In many works 
[30l - l32j the folding procedure was explored for the real 
part of the OP. Within the latter procedure the direct 
and exchange parts of the ReOP with effective nucleon- 
nucleon forces are calculated. At the same time the P 
is usually taken in a phenomenological form. Many suc¬ 
cessful applications of this model have been made for the 
proton- and nucleus-nucleus collisions (see, e.g., cycles 
of works 31, |33[ |3J|). This model was also explored in 
Refs. [35l. l36j for scattering of 10,11 Be-|-p, where the ex¬ 
change part of the folded ReOP is taken in the form of 
the zero-range prescription and the density distribution 
of 11 Be has the Gaussian-oscillator form. In the recent 
work [37] the authors account for the full exchange part 
of the ReOP, while the ImOP was calculated using the 
folding HEA formula from Refs. [38, ,39]. In Refs. [3n.l40| 
a surface term to the ImOP was added to improve the 
agreement with the data at lower energies. 

In our present work, as well as in our previous works 
considering processes with exotic He and Li isotopes jUl- 
|44| , we use microscopically calculated OPs within the hy¬ 
brid model (3^. In the latter the ReOP is calculated by 
a folding of a nuclear density and the effective NN po¬ 
tentials [32| (see also SI) and includes both direct and 
exchange parts. The ImOP is obtained within the HEA 
model [4611471 ]. There are only two or three fitting param¬ 
eters in the hybrid model that are related to the depths 
of the ReOP, ImOP and the spin-orbit part of the OP. 
Along with some phenomenological density distributions 
for He and Li isotopes we have used in our works realistic 
microscopic density obtained within the large-scale shell 
model (LSSM) [ 4 ^^]. In the present work devoted to 
processes with 10 ' n Be nuclei we use the density distri¬ 
bution for 10 Be obtained within the QMC model [s3, f5lj | 
and also the densities of 10 Be and 11 Be obtained within 
the GCM H. 

The main aim of our work is twofold. First, we 
study the elastic scattering of the neutron-rich exotic 
10 Be and n Be nuclei on protons and nuclei at energies 
E < 100 MeV/nucleon using microscopically calculated 
in our work real and imaginary parts of the optical po¬ 
tentials. Second, we estimate important characteristics 
of the reactions with 11 Be, such as the breakup cross 
sections and momentum distributions of fragments in 
breakup processes. To this end we use the model in which 
n Be consists of a core of 10 Be and a halo formed by a 
motion of a neutron in its periphery (e.g., Refs. [5314551 ]). 
The latter model is justified by the small separation en¬ 
ergy S n = 504 ± 6 KeV of a neutron from the ground 
S 1/2 state of 11 Be jUj and on the observed quite large 
total interaction cross sections of 11 Be with target nu¬ 


clei caused by the main contribution from the breakup of 
n Be on 10 Be and a neutron. The important role of the 
periphery is confirmed also by the experiments on scat¬ 
tering of 11 Be on the heavy nucleus of 208 Pb [12, where 
the prevailing mechanism is the direct breakup due to the 
long-range Coulomb force of the nucleus. Also we should 
mention the important observation of the narrow peak in 
the momentum distribution of the 10 Be fragments at the 
breakup of 11 Be scattering on the 12 C nucleus [12, that 
is, as mentioned above, a consequence of the large exten¬ 
sion of the wave function of the relative motion in the 
10 Be+n system related to the small neutron separation 
energy. By means of such a cluster model of 11 Be one can 
calculate the OPs for scattering of n Be on protons or nu¬ 
clear targets. To this end one should use the known n+p 
potential and calculate using the microscopic model the 
optical potentials of 10 Be+p (or 10 Be+A and the n+A 
potentials). Then the sum of these potentials are folded 
with a density probability of the relative motion of the 
core 10 Be and the neutron. Also, in the framework of 
this cluster model one can calculate the momentum dis¬ 
tribution of 10 Be fragments from the breakup reactions 
11 Be+ 9 Be, n Be+ 93 Nb, n Be+ 181 Ta, and n Be+ 238 U for 
which experimental data are available. 

The structure of the paper is as follows. The theoreti¬ 
cal scheme to calculate microscopically within the hybrid 
model the ReOP, ImOP, the spin-orbit part of the OP, 
the surface component of OP, as well as the results of 
the calculations of the elastic scattering cross sections of 
10,1 -^e+p and 10 ’ n Be+ 12 C are presented in Sect. [TTJ The 
basic expressions to estimate the breakup of 11 Be and to 
calculate the cross sections and the fragment momentum 
distributions of 10 Be in the diffraction and stripping pro¬ 
cesses of n Be on 9 Be, 93 Nb, 181 Ta, and 238 U are given in 
Sect, m The summary and conclusions of the work are 
included in Sect. m 

II. ELASTIC SCATTERING OF 1011 BE ON 
PROTONS AND 12 C AT E < 100 MEV/NUCLEON 

A. Hybrid model of the microscopic optical 
potential 

In the present work we calculate the microscopic OP 
that contains the volume real ( V F ) and imaginary parts 
(W), and the spin-orbit interaction ( V ls ). This OP is 
used for calculations of elastic scattering differential cross 
sections. We introduce a set of weighting coefficients Nr, 
Nj , JVjf and Nj s that are related to the depths of the cor¬ 
responding parts of the OP and are obtained by a fitting 
procedure to the available experimental data. Details of 
the constructing of the OP are given in Reis. 13(14321 . It 4 . 
The OP has the form: 


R R r dr 


. iN\ s W l I s - C ^j^- 


(f- s), 


U(r) = N R V F {r) + iN T W(r) - 2A 2 


r dr 


(1) 
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where 2A 2 = 4 fm 2 with the squared pion Compton wave 
length A 2 = 2 fm 2 . Let us denote the values of the ReOP 
and ImOP at r = 0 by V R (= V F {r = 0)) and Wi(= 
W(r = 0)). We note that the spin-orbit part of the OP 
contains real and imaginary terms with the parameters 
Vjf and W\ s related to V R and Wi by the V £ s = Vr/ 4 and 


Wj s = Wj /4. correspondingly. Here V R and Wi (and 
and W l T s ) have to be negative. The ReOP V F (r) is a sum 
of isoscalar {V r F s ) and isovector (V/y) components and 
each of them has its direct (V/g and Vfy) and exchanged 
( V F S X and V/y X ) parts. 

The isoscalar component has the form 


Vfs(r) = Vi D s(r) + Vi E s x {r) = j d 3 r p d 3 r t {p p (r p )p t (r t )v% N {s) + p p (r p ,r p + s)p t (r u r t -s)v%%(s) exp[*K(r)s/M]}, (2) 


where s = r + r t — r p is the vector between two nucleons, 
one of which belongs to the projectile and another one to 
the target nucleus. 

In the first term of the right-hand side of Eq. 0 the 
densities of the incident particle p p and the target nu¬ 
cleus pt are sums of the proton and neutron densities. 
In the second term p p and p t are the corresponding one- 
body density matrices. In our work we use for them the 
approximations for the knock-on exchange term of the 
folded potential from Refs. [58, [59] (see also jU S3). 
In Eq. 0 K(r) is the local momentum of the nucleus- 
nucleus relative motion and v F N and v FX are the di¬ 
rect and exchange effective NN potentials. They con¬ 
tain an energy dependence usually taken in the form 
g{E) = 1 — 0.003E and a density dependence with the 
form for the CDM3Y6 effective Paris potential [32j 


F(p) = C 


1 + ae~ 0p{r) 


7p(r) 


(3) 


W H in Eq. ( 0])) one can use the definition of the eikonal 
phase as an integral of the nucleon-nucleus potential over 
the trajectory of the straight-line propagation and has to 
compare it with the corresponding Glauber expression 
for the phase in the optical limit approximation. In the 
suggested scheme we use the nuclear densities and NN 
cross sections known from other sources and also the al¬ 
ready used NN potentials and amplitudes. In this way, 
the only free parameters in our approach are the param¬ 
eters Ns that renormalize the depths of the OPs com¬ 
ponents. In the spin-orbit parts of the OP the functions 
fi( r) (i = R, I) correspond to WS forms of the potentials 
with parameters of the real and imaginary parts V R , Wi , 
Ri, at [f R (r, Rr, a R ) and //(r, f?/, a/)], as they are used 
in the DWUCK4 code [6lJ and applied for numerical cal¬ 
culations. We determine the values of these parameters 
by fitting the WS potentials to the microscopically cal¬ 
culated potentials V F (r) and W(r). 


with £7=0.2658, a=3.8033, /3=1.4099 fm 3 , and 7=4.0 
fm 3 . The effective NN interactions v F N and v EX have 
their isoscalar and isovector components in the form of 
M3Y interaction obtained within o-matrix calculations 
using the Paris NN potential [3l], [32] ■ The isovector com¬ 
ponents V/y of the ReOP can be obtained by exchanging 
in Eq. @ the sum of the proton and neutron densities 
in p p (r\ by their difference and using the isovector parts 
of the effective NN interaction. In the case of the proton 
scattering on nuclei Eq. © contains only the density of 
the target nucleus. 

The ImOP can be chosen either to be in the form of 
the microscopically calculated V F (W = V F ) or in the 
form W H obtained in Ref. (38[, H2] within the HEA of the 
scattering theory [ifl [43] : 

- rp poo 

wH (r) = --JL— j 0 {kr)p p (q)p t (q)f N {q)q 2 dq.( 4) 

In Eq. 0 p(q) are the corresponding formfactors of the 
nuclear densities, /jv(q) is the amplitude of the NN scat¬ 
tering and dN is the averaged over the isospin of the 
nucleus total NN scattering cross section that depends 
on the energy. The parametrization of the latter de¬ 
pendence can be seen, e.g., in Refs. [U[6(|- We note 
that to obtain the HEA OP (with its imaginary part 


B. Results of calculations of elastic scattering cross 
sections 

In the calculations of the microscopic OPs for the scat¬ 
tering of 10 ’ n Be on protons and nuclei we used realistic 
density distributions of 10 Be calculated within the quan¬ 
tum Monte Carlo model [50, 51] and of 10,11 Be from the 
generator coordinate method [52|. In general, the QMC 
methods include both variational and Greens function 
Monte Carlo methods. In our case within the QMC 
method the proton and neutron densities of 10 Be have 
been computed with the AV18+IL7 Hamiltonian 1511 . 
As far as the GCM densities are concerned, in Ref. [52[ 
the 10 Be wave functions are defined in the harmonic os¬ 
cillator model with all p-shcll configurations. The n Be 
wave functions are described in terms of cluster wave 
functions, relative to 10 Be and to the external neutron. 
Thus, both microscopic densities effectively account for 
the non-ordinary nuclear structure peculiarities of 10,11 Be 
[26l [27j and their use is physically justified. The QMC 
and GCM densities are given in Fig. [TJ It can be seen 
that they have been calculated with enough accuracy up 
to distances much larger than the nuclear radius. In 
both methods the densities of 10 Be occur quite simi¬ 
lar up to r ~ 3.5 fm and a difference between them is 
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seen in their asymptotics. In the calculations of the OPs 
for 10,11 Be+ 12 C the density of 12 C was taken in sym¬ 
metrized Fermi form with radius and diffuseness param¬ 
eters c = 3.593 fm and a = 0.493 frn, respectively ;62j . 
The results of the calculations are compared with the 
available experimental data. All calculations of elastic 
scattering using the obtained OPs are performed by us¬ 
ing the DWUCK4 code jlili . 

1. Elastic scattering cross sections of 10,11 Be+p 

On the basis on the scheme presented in subsection 
II.A. we calculated the elastic scattering cross sections of 
10 ' n Be-|-p and compared them with the available exper¬ 
imental data. 

It is accepted that the elastic scattering of light nu¬ 
clei is rather sensitive to their periphery, where transfer 
and breakup processes also take place. Therefore, investi¬ 
gating the elastic scattering, one must bear in mind that 
virtual non-elastic contributions can also take part in the 
process. It has been pointed out in our previous papers 
|42],[43j], as well as in Refs. [M, L 3Zl, 53 , that the inclusion 
of a surface imaginary term to the OP [Eq. JTJ] leads to a 
better agreement with the experimental data. As known, 
this contribution can be considered as the so-called dy¬ 
namical polarization potential, which allows one to sim¬ 
ulate the surface effects caused by the latter. In fact, 
the imaginary part of the Is term in our OP [see Eq. ©] 
plays effectively this role. However, sometimes one needs 
to increase the absorption in the surface region and thus, 
one adds a derivative of the ImOP (surface term): 

W sf (r) = - iN s /r dW } r \ (5) 

dr 

where is also a fitting parameter. 

The results for the elastic 10 Be+p and 11 Be+p scat¬ 
tering cross sections are given in Figs. [2] and [3l re¬ 
spectively, and ^compared with the data at energies 39.1 
MeV/nucleon [63] and 59.4 MeV/nucleon j28j for 10 Be 
and 38.4 MeV/nucleon [63} and 49.3 MeV/nucleon [28| 
(see also 0 ) for n Be. In general, our analysis points out 
that more successful results are obtained in the case when 
the ImOP is taken from HE A: W{r) = W H (r) [Eq. Ipf|)]. 
We note that in the fitting procedure of the theoretical 
results to the data for elastic scattering cross sections for 
10,11 Be-|-p (and also for 10 ’ 1:L Be-t- 12 C) there arises an am¬ 
biguity in the choice of the optimal curve among many 
of them that are close to the experimental data. Due 
to this we impose a physical constraint, namely choosing 
those ReOP and ImOP that give volume integrals which 
have a correct dependence on the energy. The volume 
integrals have the forms 

ME) = J drr 2 [N R V F (r)), (6) 
J w(E) = J drr 2 [JW(r)], (7) 




47T 

ApA$ 


drr 2 


AbIF(r) - Nj f r dW ^ 
dr 


(?) 

where A p and A t are the mass numbers of the projectile 
and the target, respectively. In Eq. ([8]) we added also the 
integral over the surface term of the OP J5]). It is known 
[65} that the volume integrals (their absolute values) for 
the ReOP decrease with the increase of the energy, while 
for the ImOP they increase up to a plateau and then de¬ 
crease. The values of the Ns parameters from the fitting 
procedure and after imposing the mentioned constraint 
are given in Table [I] It can be seen that the tendency 
(the decrease of Jy and the increase of Jw) is generally 
confirmed. 


The calculated differential cross sections of 10 Be+p 
elastic scattering at energies 39.1 MeV/nucleon and 59.4 
MeV/nucleon are presented in Fig. [2] First, it is seen 
from the upper panel that the inclusion of only the vol¬ 
ume OP is not enough to reproduce reasonably well the 
data in the small angles region. Then, after adding the 
spin-orbit component to the OP the agreement with the 
data becomes better, in particular for the angular dis¬ 
tributions calculated using the GCM density at energies 
39.1 MeV/nucleon and 59.4 MeV/nucleon for angles less 
than 20° and 30°, correspondingly, as illustrated in the 
middle panel of Fig. [2] However, a discrepancy at larger 
angles remains. At the same time for the cross sections 
with the account for the Is interaction and using the 
QMC density we obtain fairly good agreement with the 
data at both energies and only a small discrepancy is 
seen at small angles at energy 59.4 MeV/nucleon. Fur¬ 
ther improvement is achieved when both Is- and surface 
terms are included in the calculations. In this case, as 
it can be seen from the bottom panel of Fig. [2j the dis¬ 
crepancy between the differential cross sections for the 
GCM density and the experimental data at larger angles 
is strongly reduced. 

In general, the account for the spin-orbit term in the 
volume OP gives a trend of an increase of the cross sec¬ 
tions at larger angles, that seems to be related with the 
change of the form of the total OP at its periphery. If 
we evaluate the quantities of the two densities of 10 Be on 
the basis of the values of the parameter N R (comparing 
which ones are closer to unity), our conclusion is that in 
the calculations without Is interaction the GCM density 
works better, while in the case with Is term in the OP 
the QMC density gives better results. A fair agreement 
between the calculated 10 Be+p angular distributions and 
the experimental data is obtained only when both Is- and 
surface contributions to the OP are included. 


In Fig. [3] are given and compared with the empirical 
data elastic cross sections for the scattering of 11 Be on 
protons at energies 38.4 and 49.3 MeV/nucleon applying 
the fitting procedure for the parameters Ns. All of them 
are calculated using GCM density of n Be. The different 
curves drawn in Fig.[5]correspond to those given in Fig. [2] 
with accounting for different contributions to the OP. 
One can see a discrepancy at small angles (9 < 30°) that 









5 



1 00 


r [fm] 



FIG. 1. (Color online) Point-proton (normalized to Z = 4) and point-neutron (normalized to TV = 6 and N = 7, respectively) 
densities of 10 Be and 11 Be obtained in the GCM [k|] and in the QMC method 5ll. PT1 . 


TABLE I. The renormalization parameters Nr, Ni, Nr, Aj s , and N^ , the total reaction cross sections ctr (in mb), and the 
volume integrals Jv, Jw\ and (in MeV.fm 3 ) as functions of the energy £=39.1 and 59.4 MeV/nucleon for the 10 Be+p 
and £=38.4 and 49.3 MeV/nucleon for the n Be+p elastic scattering 


Nucleus 

Model 

£ 

Nr 

Ni 

Nr 

n 1 / 

N 3f 

<JR 

Jv 

r(a) 

J W 

j(b ) 

J W 

10 Be 

GCM 

39.1 

0.983 

0.267 

0.000 

0.000 

0.000 

292.12 

389.408 

116.600 

116.600 

without Is and 

QMC 


1.153 

0.295 

0.000 

0.000 

0.000 

311.36 

411.344 

130.806 

130.806 

surface terms 

GCM 

59.4 

1.001 

0.802 

0.000 

0.000 

0.000 

341.18 

333.739 

263.540 

263.540 


QMC 


1.188 

0.856 

0.000 

0.000 

0.000 

356.98 

354.606 

283.464 

283.464 

10 Bc 

GCM 

39.1 

1.493 

0.492 

1.000 

0.476 

0.000 

372.50 

591.440 

216.480 

216.480 

with Is and 

QMC 


1.163 

0.318 

0.557 0.000 

0.000 

323.96 

414.911 

141.004 

141.004 

without surface 

GCM 

59.4 

1.294 

0.804 

0.190 

0.000 

0.000 

355.29 

431.427 

264.197 

264.197 

terms 

QMC 


1.014 

0.527 

0.940 

0.000 

0.000 

287.68 

302.669 

174.516 

174.516 

10 Bc 

GCM 

39.1 

0.995 

0.266 

0.095 

0.082 

0.004 

298.65 

394.161 

117.040 

122.321 

with Is and 

QMC 


1.194 

0.260 

0.075 

0.025 

0.018 

333.71 

425.971 

115.286 

139.235 

surface terms 

GCM 

59.4 

0.970 

0.000 

0.365 

1.000 

0.373 

400.26 

323.404 

0.000 

367.802 


QMC 


1.043 

0.281 

0.000 

1.000 

0.270 

389.27 

311.325 

93.053 

361.343 

11 Bc 

GCM 

38.4 

0.824 

0.659 

0.000 

0.000 

0.000 

459.05 

339.388 

293.493 

293.493 

without Is and 


49.3 

0.793 

0.805 

0.000 

0.000 

0.000 

423.52 

296.301 

301.184 

301.184 

surface terms 












n Bc 

GCM 

38.4 

0.787 

0.799 

0.000 

0.507 0.000 

458.63 

324.148 

355.844 

355.844 

with Is and 
without surface 


49.3 

0.793 

0.867 

0.123 

0.316 

0.000 

426.85 

296.301 

301.184 

301.184 

terms 












11 Be with 

GCM 

38.4 

0.849 

0.106 

0.102 

0.380 

0.152 

493.01 

349.685 

47.208 

269.903 

Is and surface 


49.3 

0.801 

0.000 

0.213 

0.394 

0.200 

436.46 

299.280 

0.000 

246.162 

terms 













seems to be related to the contributions from the surface 
region of interactions, where breakup processes play an 
important role. Similarly to the results for the 10 Be+p 
elastic scattering cross sections (see Fig. 0, the account 
for both spin-orbit and surface terms to the OP leads to 
a better agreement with the n Be+p data in the region 
of small angles. In Table Q] are given the corresponding 


values of the parameters Nr and Nj whose values deviate 
from unity of about 20 — 30% that points out that the 
hybrid model for the 0*P can be used successfully in 
such calculations. 

We would like to emphasize the fact that when con¬ 
sidering the case of the total OP [Eqs. (0) and (0], the 
values of the parameters Nj drop down sufficiently in 
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FIG. 2. 10 Be+p elastic scattering cross sections. Upper panel: without Is term; middle panel: with Is term; bottom panel: with 
both Is and surface terms. Solid lines: calculations with GCM density of 10 Be; dashed lines: calculations with QMC density of 
10 Be. Experimental data for 39.1 MeV/nucleon and 59.4 MeV/nucleon are taken from Refs. [63] and [28], respectively. 
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FIG. 3. 11 Be+p elastic scattering cross sections. Calculations are performed with GCM density of n Bc. Solid line: OP with 
both Is and surface terms [Eqs. © and ©]; dashed line: OP with Is term [Eq. ©]; dotted line: the volume part of OP from 
Eq. ©. Experimental data for 38.4 MeV/nucleon and 49.3 MeV/nucleon are taken from Refs. 63] and [28}, respectively. 




comparison with their values coming out from the two 
other cases. They are compensated in the most cases by 
the non-zero values of IVj^, iVj s , and Nj parameters. 
Here we would like to note that the Is term used in our 
calculations (with both real and imaginary parts) plays 
the similar role as the surface term applied in Ref. [371 ]. 
where, however, the imaginary Is term is disregarded. 
From our analysis made for the elastic scattering of 10 Be 
and 11 Be on protons we conclude also that the surface 
imaginary part of the OP is less necessary to fit the data 
of proton elastic scattering on the stable nucleus 10 Be, 
but it is important to give an agreement with the pro¬ 
ton elastic-scattering data of the halo nucleus 11 Be. This 
is mainly due to the specific halo structure of the 11 Be 
density distribution and its large rms radius. 

For a more complete analysis of the elastic scattering 
cross sections, we extend the incident energy region to 
lower energies on the example of the scattering of 10 Be 
on protons that has been recently studied by Schmitt et 
al. ]6fj. Moreover, this could be a test of our hybrid 
model at low energies. In Ref. [66} proton energies of 6, 
7.5, 9, and 10.7 MeV were selected to measure the elas¬ 
tic scattering cross sections for protons with 10 Be beams 
in inverse kinematics in order to provide constraints on 
optical potentials for reaction studies with light neutron- 
rich nuclei. The calculated results for the differential 
cross sections, shown as a ratio to Rutherford scattering, 
are given and compared with the data [66[ in Fig. 0] for 
energies of 7.5 and 10.7 MeV. The values of the Ns pa¬ 
rameters from the fitting procedure and the correspond¬ 
ing total reaction cross sections and volume integrals are 
listed in Table El The results shown in Fig. [I] when in¬ 
cluding in the calculations only the Is term demonstrate 
a fairly good agreement with the data. The values of 
the parameters Nr deduced from the fitting procedure 
for both energies in the case of GCM density of 10 Be 


are quite large that indicates for the specific peculiarities 
of the elastic scattering at low energies with account for 
the spin-orbit term. We also calculated the 10 Be+p elas¬ 
tic scattering cross sections at the same proton energies 
taking into account the surface term [Eq. (|5jl]. In this 
case, only the QMC density of 10 Be was tested that has 
been also used in Ref. [37], where the two other energies 
of 6 and 9 MeV were considered. The results illustrate 
that the inclusion of the surface contribution does not 
affect the good agreement obtained without it. Here we 
note that in Ref. [66| no single optical potential had been 
found to reproduce well the proton elastic scattering data 
over this range of energies. At the same time, it is seen 
from the left panel of Fig. 0] a deviation of the results of 
our model with both densities beyond 55°. Therefore, it 
would be desirable to measure the elastic channel in this 
angular range to constrain the p- 10 Be optical potential. 


2. Elastic scattering cross sections of 10,11 Be+ 12 C 

The calculated within the hybrid model elastic scat¬ 
tering cross sections of 10,11 Be-|- 12 C (their ratios to the 
Rutherford one) at the same energies, as for 10 ' n Be-|-p 
scattering, are given in Figs. [5] and [6] and compared with 
the experimental data (see also [64}). In comparison with 
the case of 10 ’ n Be-|-p, the experimental data [2^, iH} 
for the scattering on 12 C demonstrate more developed 
diffractional picture on the basis of the stronger influ¬ 
ence of the Coulomb field. It can be seen in Fig. [5] that 
in both cases of calculations of OPs with QMC or GCM 
densities the results are in a good agreement with the 
available data. It is seen also from the figures that it 
is difficult to determine the advantage of the use for the 
ImOP W = W H or W = V F , because the differences 
between the theoretical results start at angles for which 










FIG. 4. 10 Be+p elastic scattering cross sections as a ratio to Rutherford scattering at proton energies of 7.5 MeV (left panel) 
and 10.7 MeV (right panel). The solid and dashed lines show the results with QMC and GCM density of 10 Be, respectively, 
and with Is term in OP. The dotted lines show the QMC results obtained by accounting for both the Is- and surface terms in 
OP. Experimental data are taken from Ref. [66| . 


TABLE II. The renormalization parameters Nr, Ni, Nr, N \ 3 , and N^ , the total reaction cross sections gr (in mb), and the 
volume integrals Jv, Jw , an d Jw (i n MeV.fm 3 ) as functions of the proton energy E= 7.5 and 10.7 MeV for the 10 Be+p elastic 
scattering 


Nucleus 

Model 

E 

Nr 

Ni 

Nr 

N 1 / 

Nj f 

OK 

Jv 

7-0) 

rO) 

J W 

10 Be 

GCM 

7.5 

2.287 

0.473 

0.000 

0.425 

0.000 

906.19 

1215.283 

527.749 

527.749 

with Is and 

QMC 


1.244 

0.056 

0.065 

0.103 

0.000 

330.03 

603.634 

62.966 

62.966 

without surface 

GCM 

10.7 2.232 

1.129 

0.000 

0.759 

0.000 

804.23 

1144.742 

1151.009 

1151.009 

terms 

QMC 


1.915 

0.247 0.963 

0.307 0.000 

722.54 

895.179 

253.766 

253.766 

10 Be 

QMC 

7.5 

1.483 

0.000 

0.442 

0.208 

0.044 

306.28 

719.605 

0.000 

148.453 

with Is and 

QMC 

10.7 

1.354 

0.098 

0.178 

1.000 

0.193 

636.50 

632.936 

100.685 

695.676 

surface terms 













the experimental data are not available. The values of 
the parameters Nr and Ni (the depths of ReOP and 
ImOP) are given in Table IIIII From the comparison of 
these values, when GCM or QMC densities are used, one 
can see that in the case of GCM densities the values of 
the parameters are closer to unity. In this way, we may 
conclude that as in the 10 Be+p case without Is term of 
OP, the GCM density can be considered as a more real¬ 
istic one. 


III. BREAKUP REACTIONS OF 11 BE 

A. The 1(, Be+n model of n Be 

In this section we consider the characteristics of 
breakup processes of the n Be nucleus, namely diffraction 
and stripping reaction cross sections and the momentum 


distributions of the fragments. We use a simple model in 
which u Be consists of a core of 10 Be and a halo of a sin¬ 
gle neutron (see, e.g., [H). In this model the density of 
10 Be has to be given. As in Sect.[TT]we use the QMC j50| 
and GCM (52| density distributions of 10 Be. The hybrid 
model is applied to calculate the OP of the interaction 
of 10 Be with the target, as well as OP for the n+target 
interaction. In the final step of the procedure the sum 
of these potentials is folded with the respective density 
distribution corresponding to the relative motion wave 
function of the clusters in n Be. The latter is obtained by 
solving the Sclirodinger equation with the Woods-Saxon 
potential for a particle with a reduced mass of two clus¬ 
ters. The parameters of the WS potentials are obtained 
by a fitting procedure, namely, to reach the neutron sepa¬ 
ration energy S n = 504±6 KeV. They have the following 
values for 2s state in which the valence neutron in 11 Be is 
mainly bound (see Refs. 0,113): R = 2.7 fm, a = 0.52 
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FIG. 5. 10 Be+ 12 C elastic scattering cross sections. Solid lines: W = W , dashed lines: W = V F . Left panel: calculations 
with GCM density of 10 Be; right panel: calculations with QMC density of 10 Be. Experimental data for 39.1 MeV/nucleon and 
59.4 MeV/nucleon are taken from Refs. [63| and [28j, respectively. 


TABLE III. The renormalization parameters Nr and Ni, the total reaction cross sections or (in mb), and the volume integrals 
Jy and Jw (in MeV.fm 3 ) as functions of the energy 15=39.1 and 59.4 MeV/nucleon for the 10 Be+ 12 C and 15=38.4 and 49.3 
MeV/nucleon for the n Be+ 12 C elastic scattering 


Nucleus 

Model 

E 

W 

Nr 

Ah 


Jv 

j( a ) 

J w 

10 Be 

GCM 

39.1 

W H 

0.939 

0.708 

104.539 

255.156 

283.037 




V F 

0.816 

0.465 

105.958 

221.733 

126.355 



59.4 

W H 

1.013 

1.010 

101.052 

238.122 

302.581 




V F 

0.884 

0.577 

102.635 

207.798 

135.633 

10 Be 

QMC 

39.1 

W H 

0.888 

0.620 

105.332 

245.613 

249.769 




V F 

0.782 

0.434 

106.878 

216.294 

120.041 



59.4 

W H 

0.970 

0.887 

101.616 

231.953 

267.782 




V F 

0.849 

0.534 

103.035 

203.019 

127.694 

11 Be 

GCM 

38.4 

W H 

0.769 

0.711 

127.123 

216.879 

287.235 




V F 

0.708 

0.521 

126.825 

199.676 

146.937 



49.3 

W H 

0.820 

0.883 

124.406 

213.754 

300.193 




V F 

0.743 

0.574 

123.302 

193.682 

149.628 


fm and Vq = 61 MeV. The rms radius of the cluster for¬ 
mation is obtained to be 6.87 fm. 

The s-state (l = 0,n = 1,2) of the relative motion of 
two clusters has the form: 


4>oq( s ) = 0o l) ( s )^==’ n= 1,2. (9) 


The corresponding density distribution is the probability 
of both clusters to be at a mutual distance s: 

pSAs)H^A)l 2 = T|4 n| Ml 2 . (io) 

Within the 10 Be+n cluster model, in order to calculate 
the 11 Be breakup in its collision with the protons and nu¬ 
clear targets, one should calculate two OPs of 10 Be+p(or 
A) and n+p(or A) scattering: 
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FIG. 6. 11 Be+ 12 C elastic scattering cross sections. Solid lines: 
W = W H ; dashed lines: W = V F . For n Be GCM density 
was used. Experimental data for 38.4 MeV/nucleon and 49.3 
MeV/nucleon are taken from Refs. [63] and fH], respectively. 


U^ b ’ n \r) = F (b ’ n) +iW {b ’ n) = JdstfHs) C/ c (n) (r + (l/ll)s) + C/^(r-(10/ll)s) =2 np^\s)s 2 ds 

/•l 

x J dx t/ c (n) (yr 2 + (ls/ 11) 2 + r( 2 /ll)sa;) + U^ ] (Vr 2 + ( 10 s/ll ) 2 - r( 20 /ll)sa;) . ( 11 ) 


In Eq. (fill) r — (10/ll)s = r„ and r + (l/ll)s = r c give 
the distances between the centers of each of the clusters 
and the target, and s = Sj + s 2 = ( 10 /ll)s + (l/ll)s 
determines the relative distance between the centers of 
the two clusters, si and S 2 are the distances between the 
centers of n Be and each of the clusters, correspondingly. 
The respective OPs for the 10 Be+A and n+A scattering 
are calculated within the microscopic model of OP from 
Sect. II.A. 

In the case of the 11 Be breakup on the proton target 
the n+p potential is taken in the form [ 68 ] (in MeV): 

U^ n) = Vn P = u(r)(! + ij). ( 12 ) 


with 


the expressions of the 5-matrix (as a function of the im¬ 
pact parameter 6 ): 


5(6) = exp 


i 

hv 



U{^b 2 + z 2 )dz , 


where 


(14) 


U=V + iW (15) 

is the OP. For negative V and W one can write 

S(b) = 


x exp 


i r°° 

/ \W\dz 

hv oo 


( 16 ) 


v(r) = 120e _1 ’ 487r2 —53.4e _0 ' 639r2 —27.55e -0 ' 4657 ' 2 , (13) 
where 7 = 0.4. 

For calculations of breakup cross sections and momen¬ 
tum distributions of fragments in the 10 Be+n breakup 
model we give here briefly the eikonal formalism, namely 


and, correspondingly, 

|5(6)| = exp 


1 r°° 

-- / \W\dz 
. hv j 


(17) 


In our case W is the imaginary part of the microscopic 
OP [Eq. (fill) ]. 1 5(6 )| 2 gives the probability that after the 
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collision with a proton (2 —> oo) (in the * 11 Be+p scatter¬ 
ing) , the cluster c or the neutron with impact parameter 
b remains in the elastic channel (i = c, n): 


probability function of the k-momentum distribution of 
a cluster in the system of two clusters as a function of 
the impact parameter b |16j: 


|5' i (6)| 2 = exp 


2 

hv 



dz 


Wi{Vb 2 + z 2 ) 


(18) 


The probability a cluster to be removed from the elas¬ 
tic channel is (1 — |5j 2 ). The probability of the case 
when both clusters (c and n) leave the elastic channel 
is (1 — |5'„| 2 )(1 — |S'c| 2 )- As shown in the next subsec¬ 
tion, Eqs. (fl4ll - (fT8ll take part in the calculations of the 
diffraction breakup and stripping reaction cross sections. 


B. Momentum distributions of fragments 

The necessary quantity to calculate the diffraction 
breakup and absorption scattering cross sections (dif¬ 
ferential and total) and momentum distributions is the 

I 


d 3 Pn( b,k) 
dk. 


(27r) 3 


ds^(s)n(b,r J _)^” ) (s) 


■ (19) 


In Eq. m ^(b, rj_) is given by the products of two S- 
functions S c and S n [Eqs. (fTUl) - ([TUl) ] of the core 10 Be and 
the neutron, ^k(s) is the continuum wave function, k is 
the relative momentum of both clusters in their center- 
of-mass frame, and the vector ri is the projection of the 
relative coordinate s between the centers of the two clus¬ 
ters on the plane normal to the 2 -axis. The bound-state 
wave function </>oo of the relative motion of two clusters is 
given for the s-state by Eq. ©. As to the wave function 
in the final state <j) k, we will neglect its distortion and, 
thus, replace it by jo(ks) in the case of the s-state. Then, 
following Ref. |16J, the probability function has the form 


d 2 Pp(b, k) 
dk^dk± 


k± 

167t 3 A; 2 



d(cos 9 a ) g(s) sin ( ks ) 


2 


c%fl(b, rj_) , 


( 20 ) 


where 


fi(b,r x ) = S c (b c )S„(b„) 


( 21 ) 


and g(s) = s0q"H s ), being given by Eq. ©. 

Hence, the diffraction breakup cross section has the 
form 


C 2 7T 


bn.dbn 




dk^dk± 


(—) = 

V dkL J diff 

( 22 ) 

In Eq. BUI) d 2 Pn(b, k)/dkLdk± is given by Eq. (BUI) . The 
integrals over b n and ip n mean integration over the impact 
parameter b„ of the neutron with respect to the target. 

The cross sections of the stripping reaction when the 
neutron leaves the elastic channel is [1(|: 

(dfc^) str = 2^i b n db n dv n [l-\S n {b n )\ 2 ] 


x / pdpdip p \S c (b c )Y 


dzcos(kLz)cj) 0 p 2 + 2 2 ^ 


.(23) 


Equation (1231) is obtained when the incident nucleus has 
spin equal to zero and for the s-state of the relative mo¬ 
tion of two clusters in the nucleus with s = r c — r n , 
p = b c — b„, s = p + z and 

b c = yj s 2 sin 2 9 + & 2 + 2 sb n sin 6 cos(<^ — ip n ) (24) 


coming from b c = h n + b, where b = s sin 9 is the pro¬ 
jection of s on the plane normal to the 2 -axis along the 
straight-line trajectory of the incident nucleus. 

In the end of this subsection we note that the real and 
imaginary parts of the OPs taking part in Eq. m and in 
the S'-matrices [Eqs. (fl4l) - (fl8l) [ are used for calculations of 
the cross sections [Eqs. (fl9l) - (BUl) [ in the cases of scatter¬ 
ing and breakup of 11 Be on protons and nuclei that will 
be considered in the following part of our work. They are 
calculated microscopically within the hybrid model given 
in subsection II.A. 


C. Results of calculations of breakup reactions 

In this subsection we perform calculations of the 
breakup cross sections of 11 Be on the target nucleus 
9 Be and heavy nuclei, such as 93 Nb, 181 Ta, and 238 U, 
and compare our results with the available experimental 
data [lJl. The densities of these heavy nuclei needed 
to compute the OPs are taken from Ref. [iiol . The 
diffraction and stripping cross sections (when a neu¬ 
tron leaves the elastic channel) for reactions n Be+ 9 Be, 

11 Be+ 93 Nb, 11 Be+ 181 Ta, and n Be+ 238 U are calculated 
from Eqs. (BUI) and (B3l) . The obtained results are illus¬ 
trated in Figs. [7J [HI [9] and [TUI respectively. We note the 
good agreement with the experimental data from light 
and heavy breakup targets. The obtained cross sections 
for the diffraction and stripping have a similar shape. 
The values of the widths are around 50 MeV in agreement 
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FIG. 7. Cross sections of diffraction breakup and stripping 
reaction in 11 Be+ 9 Be scattering at E = 63 MeV/nucleon. 
Experimental data are taken from Ref. EH- 



k L [MeV/c] 


FIG. 8. The same as Fig. 0 but for 11 Be+ 93 Nb scattering. 


with the experimental ones. Our results confirm the ob¬ 
servations (e.g., in Refs. 00 ) that the width almost 
does not depend on the mass of the target and as a re¬ 
sult, it gives information basically about the momentum 
distributions of two clusters. Here we note that due to 
the arbitrary units of the measured cross sections of the 
considered processes it was not necessary to renormal¬ 
ize the depths of our OPs of the fragments-target nuclei 
interactions. 


IV. CONCLUSIONS 

In the present work the hybrid model is applied to 
study characteristics of the processes of scattering and 
reactions of 10 Be and n Be on protons and nuclei. In 
the model, the ReOP is calculated microscopically in a 
folding procedure of the densities of the projectile and 
the target with effective NN interactions related to the 
^-matrix obtained on the basis of the Paris NN potential. 
The ReOP includes both the direct and exchange terms. 
The ImOP is calculated microscopically as the folding OP 
that reproduces the phase of scattering obtained in the 
high-energy approximation. The only free parameters in 



FIG. 9. The same as Fig. [7] but for 11 Be+ 181 Ta scattering. 



FIG. 10. The same as Fig. 0 but for 11 Be+ 238 U scattering. 


the hybrid model (Ns) are the coefficients that correct 
the depths of the ReOP, ImOP, and the spin-orbit parts 
of OP. Their values are obtained by a fitting procedure to 
the experimental data whenever they exist. Additionally, 
in some cases the surface absorption is accounted for by 
including another term to the OP that requires one more 
fitting parameter. The density distributions of 10 Be ob¬ 
tained within GCM and QMC microscopic methods and 
of n Be from GCM are used. The resulting within the hy¬ 
brid model OPs are applied to calculate characteristics 
of various processes. 

The results of the present work can be summarized as 
follows. 

(i) Elastic scattering cross sections of 10 Be and n Be 
on protons and 12 C are calculated using the microscopic 
OPs for energies E < 100 MeV/nucleon and compared 
with the existing experimental data. In order to resolve 
the ambiguities of the magnitudes of the depths of the 
OPs, the well established energy dependence of the re¬ 
spective volume integrals of the OPs is taken into ac¬ 
count. The theoretical approach gives a good explana¬ 
tion of a wide range of empirical data on the 10 ’ n Be+p 
and 10 ’ n Be+ 12 C elastic scattering. It was established 
that the obtained by fitting procedure values of the co¬ 
efficients Nr (depths of ReOP) are close to unity. The 
correction of the ImOP by factor Nj is in some cases 
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larger, e.g., for 1 2 3 4 5 6 7 8 9 10 Be+j» at energy E = 39.1 MeV/nucleon 
in the case when the spin-orbit (Is) component is not ac¬ 
counted for. The inclusion of a surface term to the OP 
leads to a better agreement with the experimental elastic 
scattering cross section data. We conclude that, in gen¬ 
eral, the hybrid model for microscopic calculations of the 
OPs gives the basic important features of the scattering 
cross sections and can be recommended and applied to 
calculate more complex processes such as breakup reac¬ 
tions, momentum distributions of fragments and others. 

(ii) Apart from the usual folding model, we use another 
folding approach to consider the 11 Be breakup by means 
of the simple 10 Be+n cluster model for the structure of 
11 Be. Within this folding model we construct the OP 
of the interaction of 10 Be with the target, as well as the 
n+target interaction. Using the cluster OPs 10 Be+p(or 
A) and n+p (or A) the corresponding functions S c and S n 
(S'-matrices) for the core and neutron within the eikonal 
formalism are obtained. 

(iii) The calculated S c and S n functions are used to get 
results for the longitudinal momentum distributions of 
10 Be fragments produced in the breakup of 11 Be on differ¬ 
ent targets. This includes the breakup reactions of 11 Be 
on 9 Be, 93 Nb, 181 Ta and 238 U at E = 63 MeV/nucleon 
for which a good agreement of our calculations for the 
diffraction and stripping reaction cross sections with the 


available experimental data exist. The obtained widths 
of about 0 MeV/c are close to the empirical ones. Fu¬ 
ture measurements of such reactions are highly desirable 
for the studies of the exotic 11 Be structure. The accu¬ 
rate interpretation of the expected data requires more re¬ 
fined theoretical approaches, for instance that of Ref. [70] 
within the CDCC method and its extensions to study the 
effects of the dynamic core excitation, especially its large 
contribution to nuclear breakup in the scattering of halo 
nuclei. 
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